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ABSTRACT 

In this talk I discuss both the present status and some recent work on the 
Kazakov-Migdal Model which was originally proposed as a soluble, large N realization of 
QCD. After a brief description of the model and a discussion of its solubility in the large 
limit I discuss several of the serious problems with the model which lead to the conclusion 
that it does not induce QCD. The model is nonetheless a very interesting example of a 
Gauge Theory and it is related to some very interesting Matrix Models. I then outline a 

technique which uses "Loop Equations" for solving such models. A Penner-like model is 
then discussed with two logarithmic singularities. This model is distinguished by the fact 
that it is exactly and explicitly soluble in spite of the fact that it is not Gaussian. It is 
shown how to analyze this model using both a technical approach and from a more physical 
point of view. 



1. Introduction 



The Kazakov-Migdal Model is an SU(N) Lattice Gauge model which was 
proposed several years ago in the hope that it would be both soluble in the large 
limit and that it would give the correct long distance behavior of QCD in the 
continuum limit. This model was thus also known as "Induced QCD". Although it is 
now widely believed that this model does not induce QCD much has been learned from 
the study of this model, both about ordinary and about gauged Matrix Models. Work 
is still in progress on a class of "soluble" , non-Gaussian Kazakov-Migdal Models. 

This talk will begin with a brief review of the Kazakov-Migdal Model (the 
KM Model) and why it was originally though that it might induce QCD. This will 
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be followed by a discussion of the solubility (i.e. tractability) of the model at large 
N. Some of the main problems with the model will then be described including a 
discussion of the extra local Zjv symmetry and the absense of critical behaviour in 
the Gaussian Model. This will then explain why it is extremely unlikely that this 
model could induce QCD. This will be followed by an outline of a powerful method 
for solving KM Models. I shall conclude this talk with a discussion of a soluble, 
non-Gaussian, ("Penner") model which is explicitly soluble in the large N limit and 
which is the focus of continued research. 

2. The Model 

The KM Model is a d-dimensional lattice Gauge Theory containing an 
SU(N) Gauge Field Uxy defined on the links of the lattice and a Hermetian (adjoint) 
scalar "Higgs" Field defined on the sites of the lattice. The model is defined in 
terms of the Euclidean Partition Function 

links sites 

with 

5 = TV ^ tr (y($,)) - TV ^ tr ($,t/,,$,C/-i) (2) 

sites,x links,xy 

where V{^) is, at this stage, an unspecified potential for the field $ and where the 
second term is the standard gauge invariant kinetic term for the adjoint field Note 
that this action is simply the lattice action for an adjoint scalar field except that the 
Wilson term which is the kinetic term for the Gauge Field is missing. This Wilson 
term 

^ Tr {UiU2U3Ua) + h.c. (3) 

plaquettes, 1234 

is missing in order to assure the solubility of the model at large N. 

The original hope was that the integration over $ would induce a kinetic 
term for the Gauge Field which would lead to a theory which, in the continuum limit, 
would be QCD (or, more precisely. Quantum Gluodynamics, since quarks have not 
been introduced at this stage). The arguments supporting this claim can be found in 
the original papers 
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3. Large N Solubility 



The missing Wilson term (3) in Eq. (2) allows for exact integration over 
all the Gauge Fields Uxy on the lattice. To see this note that each term in the action 
contains one and only one Uxy and each such Gauge Field appears in only one term 
in the action. The integral over U can be done by using the formula'*' 

where 

A(</>) = n ('^^ - (5) 

i<j 

is the Vandermonde determinant for 0. This formula can be applied to the integral 
in Eq. (1) by diagonalizing the Hermetian Matrices and which appear in Eq. 
(2). 

The main point is that the result of the integral depends only on the 
eigenvalues of the $'s. The integral over each matrix $ (which is an integral with the 
Hermetian measure) can now be written as an integral over the eigenvalues (pi of $ 
since 

Jv^-- = j l[d(l>iA\(l>) ■ ■■ (6) 

i 

when the integrand depends only on the eigenvalues of $. It follows that the Partition 
Function, Eq. (1), is of the form 

Z = J l[d(j)x,i ^\(px) exp[-,Seff] (7) 

where S'eff depends only on the cigcnvahies of the $'s. In fact it contains one term 
which involves the (f)x,i at neighbouring sites and a potential term involving only one 
site at a time. 

In the limit of large the above integral can be done by finding the Clas- 
sical Minimum of the action S^s- This is assumed to occur when the eigenvalues of 
^x a-re independent of x. The basic problem is then to minimize the effective action 



(8) 



with respect to all the 0^. The result will be a specific set of eigenvalues 0i...0jv which, 
in the limit of large N, is described by a density of eigenvalues (proportional to 
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the number of eigenvalues in the vicinity of 0) which is normalized so that J dXp{X) = 
1. This minimization can be carried out in practice for only a very limited number 
of special choices for the potential V. We shall discuss this further later in this talk. 



4. Problems with "Induced QCD" 

Despite the solubihty of the KM Model at large N, several serious problems 
were recognized immediately with the hypothesis that the KM Model induces QCD. 
The first of these problems is the presence of an additional local Zj^ symmetry in the 
model . The absense of the Wilson term allows us to multiply the Gauge Field Uxy 
on any link (xy) by an element of the center of the Gauge Group without affecting 
the action. More formally, the transformation 

Uxy exp (^^^ Uxy (9) 

(which is also an element of SU(N) if k is an integer) leaves the action Eq. (2) 
invariant even if a different k is chosen for each link (xy) . (A Wilson term would ruin 
this invariance.) 

The main problem with this symmetry is that it implies the vanishing of 
all Wilson Loops since 

<Uu>=^-v(^)<]lu> (10) 

loop loop 



if we multiply any one of the links in the loop by exp (^). This forces all Wilson 

Loops to vanish. Furthermore, since the symmetry is local, it cannot be broken'^' . 
This is a disaster if the model is to induce QCD since, in QCD, we expect an area 
law for Wilson Loops. 

Several proposals have been made to solve this problem, none of which 
have been particularly successful. Attempts at solving this problem have included 
the introduction of Fundamental Representation Fermions or a subleading (in 1/A^) 

Wilson term both of which break the local Zn symmetry at the expense of the exact 
solubility of the model. It has also been suggested that "Filled Wilson Loops" which 
are invariant under the local Z^ replace the ordinary Wilson Loops in this model. 

A second very serious problem with the idea of inducing QCD from the 
KM Model was discovered when Gross solved the Gaussian KM Model, i.e. the 
model with V{(f)) — m^(fP/2, explicitly. The argument is as follows. It turns out 
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that the perturbative mass of the "Higgs" scalar is — 2d {d is the dimension of 
the spacetime.) The most convincing arguments that the KM Model induced QCD 
were made for precisely this Gaussian potential for which it was expected that the 
model would have critical behaviour and a continuum limit when 2d,. The 

exact solution to the model proved that no such critical behaviour was present and 
thus the Gaussian model certainly did not induce QCD. 

This, by itself, does not rule out that some non-Gaussian model may induce 
QCD but, in light of the fact that the main argument that the model induced QCD 
did work for the Gaussian case and in light of the Zj^ problem it seems very unlikely 
that any non-Gaussian model would induce QCD. There are, in fact, several other 
problems with the idea of Induced QCD which are discussed in the literature. 

Despite the failure of the model to induce QCD, the KM Model is still a 
very interesting Gauge Theory and it is very interesting to find non-Gaussian models 
which can be solved exactly. Fortunately there is a Penner-hke model with a loga- 
rithmic potential which can be solved explicitly (Ref. 1) and which will be discussed 
below. 



5. A Technique for Solving KM Models 



In this section I shall review briefly the method of Dobroliubov, Makeenko 
and Semenoff (Ref. 1) for finding the extremum of the action (8). The basic idea is 
to begin with the original action (2) and to consider the following two quantities: 



(11) 



and 



G(A,^) = (^Tr 



■ — U~' 



(12) 



where the average is with respect to the original Path Integral and where the expres- 
sion of ^^(A) in terms of the density of eigenvalues p is valid in the large N limit. 

The main idea is to write a set of equations, similar to Schwinger-Dyson 
Equations and often called "Loop Equations" for these quantities. These equations 
are of the form 



VUV^ 



d 



d^ 



X,IJ 



1 



A-$. 







and 



/ 



VUV^ 



d^ 



X,IJ 



A-$. 



-a 



xy- 



'- U~' 



(13) 



(14) 
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The first equation Eq. (13) will yield an equation for E{X) and G(A, z/) whereas the 
second equation will, in general, involve quantities with more factors of 1/(A — 0). 

In solving these equations one needs to make heavy use of the large N limit 
in which the integral is dominated entirely by a single matrix $ which minimizes the 
effective action Eq. (8). Even averages over the Gauge fields U which appear in Eq. 
(14) are performed in this background $ field. In fact it is useful to treat 

as an unknown whose only nonvanishing components are defined as 

Q^ = ^ y D[/|[/,,f (16) 

where $ is the minimum of the effective action which has yet to be determined. 
Obviously Cy depends on It is also useful to define a quantity 

Ai = Cy0j (17) 

where (pj are the eigenvalues of Now the unknown but fixed values of allow us 
to identify an integer with every eigenvalue of $. We can thus replace the quantity 
Aj with a quantity A((/)) which is defined via the correspondence of with i. 

The quantities defined above are very useful in solving the "loop equations" 
It turns out that there arc several cases in which these loop equations close and one 
can solve explicitly for -E'(A) and G(A, z/). In these cases the density of eigenvalues 
can be extracted from E{X) using Eq. (11) which imphes that E{X) has a branch 
cut along the support of p whose discontinuity is proportional to p. Cij can also be 
determined as the double discontinuity across the cut of G{X, u) in A and in //. 

The simplest case in which these equations close is the Gaussian case in 
which V{(t)) = m202/2. This case is discussed in detail in Ref. 1. What is even more 
interesting is that this system can be solved for a non-Gaussian, though singular 
potential of the Penner'^"' (logarithmic) type. An outhne of the solution follows. 

Begin with the "Loop Equations" Eqs. (13) and (14). The first equation 
Eq. (13) leads to the equation 

E(Af-(m:^^) =0 (18) 

Here V' — dV/dcf) The second equation Eq. (14) can be simplified if we assume a 
simple form for the quantity V' — 2{d — 1)A. The soluble "Penner" case consists of 
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requiring this quantity to have a simple pole at some point ^ so that 

V'{<t>) - 2{d - l)m + B (19) 

This "ansatz" allows a significant simplification of the equations. Using the asymp- 
totic conditions 

E{X)r.^ + ^ + ... (20) 

G{X,u)^^ + ... (21) 

one eventually derives a quadratic equation for E{X) in terms of a single, unknown 
quantity 4> which must be determined self consistently. This quadratic equation for 
E{X) is simply 

E{Xf [(A -B){X- 0] + ^(A) [(A -B){X- Oii - B) + q{B - - (A - 0] 

+ [{X + ^){B-O + {e-B^)]=0 (22) 

Besides allowing us to solve for E(X) and thus for the density of eigenvalues 
p(A) we can now find the potential which leads to the ansatz (19). It is not difficult 
to show that if we were to substitute the expression 

V'{X) - 2dA{X) = + + {^-B) (23) 

into equation (18) we would obtain precisely the same quadratic equation (Eq. (22)) 
for E{X) as Equation (22). It thus follows that 

A(A) = + e (24) 

and thus the potential V is given by 

^,(,)^_^^(2rf^iK|zi)+B+p,_i)5 (25) 

We are thus able to solve the KM Model for a potential with two loga- 
rithmic singularities plus a linear term when the various coefficients are related as 



1 There is actually a subtle assumption being made here. It is not guaranteed that V will simply 
be the integral of V' unless V is of the special form Tr/(0) where / is some ordinary function 
of (p. The point is that since we have made an ansatz for V' — {2d — 1)A rather than for V we 
are not guaranteed that V will be of this simple form. 
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in Eq. (25). 



6. Relationship with an Ordinary Matrix Model 

We are now ready to sec that the solution to the above KM Model is closely 
related to the solution of an ordinary Matrix model with the potential V' — 2dA given 
by Eq. (23). To see this consider an ordinary matrix model defined as an integral 
over a single N x N Hermetian matrix $ 

Z = J P$exp [-NTrW{^)] (26) 

This model can be solved by defining -E'(A) in analogy with Eq. (11) as 



E(X) = < ^ > (27) 



and writing a "loop equation" in analogy with Eq. (13)which leads to an equation 
identical to Eq. (18) 

,2 /W'W 



It follows that if, for W', we choose V' — 2dA given by Eq. (23) we arrive at precisely 
the same equation for E{X) and thus the same distribution of eigenvalues p(A). 

We are thus led to consider an ordinary Matrix Model with two logarithmic 
singularities so that W is given by 

W\\) = — ^ + + C (29) 

A - 771 A - ?72 

This model is related to the KM Model above provided (see Eq. (23)) 

ri + r2 = 1; C = 772 - ?7i (30) 



It is most interesting to consider the general two pole Penner Model above 
and to discuss the KM case as the special case when Eq. (30) is satisfied. In the 
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general case Eq. (28) can be written as: 



E{Xf - W'{X)E{X) + J2 = (31) 

i=l 

Using the asymptotic expansion Eq. (20) one can trade the two unknowns E{r]i) and 
E{r]2) for (j) via the relations 

2 

Y,nE{r^i) = C 
1=1 



^ riE{r]i)r]i = n + rs - 1 + C0 (32) 
1=1 



The solution for E{X) is given by 



E{X) = \ I W'{X) ± 



X^W- 41:^^1 (33) 
\ i=i ^ 



Now recall that E{X) = J d(f)p{(f))/{X — 0) has a branch cut precisely on 
the support of p. Furthermore on such a branch cut of E{X) 

E(X + ie) - E(X - ie) = 27rip(A) (34) 

Thus the first step in finding the density of eigenvalues is to find the location of the 
branch points of E{X) which can be done by finding the places where the term inside 
the square root vanishes. This is a quartic equation and thus, in general, there will 
be two branch cuts. So the general solution will be a two cut solution. Of course 
there will be cases when there are no real solutions to the quartic equation and there 
will be cases when wc have a one cut solution. 

The technical details of the various solutions for various values of the pa- 
rameters will be discussed in more detailed papers on this subject.'"' What I would 
like to discuss in this talk is a general physical picture of what kind of solutions one 
should expect. Before doing so, however, there is one further important technical 
point. When any of the rj are positive the potential has an infinitely attractive log- 

[12] 

arithmic well. In such cases it is known from the single pole Pcnncr Model that 
the branch cuts may have to circle the singularity. Mathematically this means that 
there is no real solution for p though physically it may be possible to interpret this 
as a "condensation of eigenvalues" at the singularity. 
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7. Physical Picture of Solutions 



It is instructive and not difficult to get a good physical picture of the 
location of the eigenvalues of $ for the solution of both the simple Matrix Model 
and the KM Model. In particular it is not difficult to decide whether solutions exist, 
whether they are unique, where the eigenvalues are relative to the singularities and 
whether there are one-cut or two-cut solutions. 

We begin with the simple Matrix Model given by Equation (26) which can 
be written, using Eq. (6) as 



(35) 



In large N, we must minimize the effective Action 



- J2 log (0z -cf>jf + Nj2 (36) 

Notice that the minimization of this Action corresponds to an analogue classical 
mechanical problem in which a large number N of particles are located at locations 
0i...0jV along a line. They are each subjected to a central potential W{(f)) and to a 
two-body repulsive potential log(0j — Thus, for example, if the potential has 

a single minimum at some point (f)o, one expects the eigenvalues to collect near this 
minimum with some distribution which is controlled by the logarithmic repulsion. If 
the potential has several minima then, presumably, there are a whole class of classical 
solutions in which a varying proportion of particles are located in each of the wells. 
In the Penner-like model in which there are at most two such minima, these various 
cases can be distinguished by the value of (p. 

Let us now apply these ideas to the potential (29). The simplest case 
occurs when both ri and r2 are negative. Recall that Eq. (29) gives the derivative of 
the potential so that in this case the potential has two repulsive singularities at rji 
and 1]2 and a linear term with either a positive or a negative coefficient. In either case 
the potential has precisely two minima. If, without loss of generality, we consider the 
case in which C>0 then one minimum is between rji and 772 smd the other is at a point 
greater than 772 (assuming 772 > ^i)- The potential is unbounded below for (f) < ri\. 
The fact that the potential is unbounded does not affect the existence of solutions. In 
fact classically the analogue "particles" can exist in the two minima without being 
significantly affected by what occurs for < 771. Thus we expect a whole class of 
solutions which can be parameterized by which are, in general, two-cut solutions, 
in which some of the eigenvalues live in a compact region between 771 and 772 and the 
remainder live in another compact region near the minimum beyond 772. 
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The next case to consider is that in which both ri and r2 are positive. In 
this case the potential has two attractive singularities at rji and at 772- It has two 
maxima but no minima. In this case there are clearly no nonsingular solutions to the 
analogue problem though, as discussed previously, there may be some solutions in 
which the cuts of E{X) circle the singularities and which can be interpreted as a con- 
densation of eigenvalues. These turn out to be cases in which the quartic form under 
the square root in Eq. (33) has real solutions but, in order to get both the asymptotic 
(large A) behaviour of E{\) and its behaviour near the singularities correctly the cuts 
must circle one or more singularities. 

The final case, in which the have opposite sign, is the most interesting 
case for the KM Model. If, without loss of generality, wc assume that r]i < r]2 and 
C > then we identify two cases. If ri < and r2 > then the potential is 
repulsive at ri and attractive at r2. In this case we expect no real solutions, though 
"condensation of eigenvalues" may occur at r2. In the case ri > and r2 < the 
potential is attractive at ri and repulsive at r2- Thus there is a minimum for > 772 
and we should expect at least one real one-cut solution to this problem for some 
specific value of 4>. 

To see that this last case is the one relevant to the KM Model recall Equa- 
tion (30) which states that for the KM case ri + r2 = 1. Now in order to get a real 
solution we need at least one of the ri to be negative. In this case we see that the 
other one must be positive so that the sum equals 1. We see from the preceding 
paragraph (and from the fact that C — 7)2 — rji > 0) that if r2 > we expect no real 
solutions whereas if r2 < we expect a single cut solution, where the cut is near the 
minimum at > 772- 

Let us now compare this with Eq. (23). We should identify 772 as rj, rji as 
B and r2 as g'. Thus we expect a real single cut solution only if < 0. The form of 
the KM potential can then be deduced from Eq. (25). Note that the qualitative form 
of V is not uniquely determined. Even though i? — ^ < 0, the sign oi B + {2d — 2)^ 
depends on the precise values of B and ^. What is however uniquely determined is 
the sign of the singularities (provided d > 1). If g < then {2d — l)(g — 1) < 
and both singularities in V as opposed to W are repulsive. The integral over the 
Gauge Fields force a renormalization of the analogue potential from V to the effective 
Matrix Model form W. 

In the previous paragraphs we have shown how to estimate the approximate 
location of the eigenvalues for an ordinary Matrix Model or for a KM Model given 
the relationship between a KM Model and an ordinary Matrix Model. It is, however 
possible to get some information from the KM Potential directly. To see this suppose 
we have a KM Model with an arbitrary potential V. Our job is to find the extrema 
of the Effective Postential given by Eq. (8). This can be rewritten in the following 
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form" 



Seff = iV^ {V{cpi) - d + (d - 1) log - <\>j)f - d log (detye^('^^-^^)'/2) 

(37) 

Using the fact that the Integral (4) is finite and nonzero when (j) = ip and when any 
two eigenvalues are equal, it is easy to see that minimization of Eq. (37) corresponds 
to the following analogue mechanical problem. We imagine again N particles living on 
a line with a central potential which isV — d(f)^ (not simply V), with an interparticle 
potential which is no longer 2-body but which has the following property: Whenever 
any two eigenvalues approach each other there is a repulsion (the Action diverges) 
and whenever any two eigenvalues are separated from each other by a large distance, 
there is an attraction which tries to bring them closer together. This is quite different 

from the physical picture which emerges from the simple Matrix Mo del 
8. Summary 

At this point there is no evidence that the Kazakov-Migdal Model pro- 
vides a soluble, large N realization of QCD and such a hypothesis has many serious 
difficulties. We have, however, learned much from the study of the KM Model both 
about Matrix Models and about this new, interesting class of Gauge Theories. 



t Unfortunately the compatibility between this picture and the one which emerges from writing 
an effective one Matrix Model and using the Physical Picture in that case, is not at all clear. 
In fact, as mentioned briefly earlier in this article, there is some possibility that the solution 
which corresponds to the one Matrix Model is not the KM Model with the simple potential 
TrV given as the integral of V but one with a more complicated matrix structure. These issues 
will be discussed in a more detailed publication. 



12 



REFERENCES 



1. DMS paper M.I. Dobroliubov, Yu. Makeenko, G.W. Semenoff, Mod. Phys. Lett. 
A8 P.2387 1993; Yu Makeenko, Mod. Phys. Lett A8 (209) 1993; Yu Makeenko, 
Mod. Phys. Lett. B422 (237) 1994 

2. V.A. Kazakov and A.A. Migdal, NucL Phys. B397 (214) 1993 

3. V.A. Kazakov and A.A. Migdal, Ref V.A. Kazakov and A.A. Migdal, Nucl. 
Phys. B397 (214) 1993; A.A. Migdal, Mod. Phys. Lett. A8 (359) 1993; A.A. 
Migdal, Mod Phys. Lett. A8 (153) 1993 

4. Harish-Chandra, Amer. J. Math. 79 (87) 1957; C.Itzykson and J.B. Zuber, J. 
Math. Phys. 21 (411) 1980; M.L. Mehta, Comm. Math. Phys. 79 (327) 1981 

5. I. Kogan, G. Semenoff, N. Weiss, Phys. Rev. Lett 69 (3435) 1992 

6. S. Elitzur, Phys. Rev. D12, 3978. 

7. A.A. Migdal, Mod. Phys. Lett. A8 (359) 1993 

8. Phys.Lett. B302 (283) 1993. 

9. D. Gross, Phys.Lett. B293 (181) 1992 

10. Yu. Makeenko, "Critical Scaling and Continuum Limits in the D^l Kazakov- 
Migdal Model" |hep-th 9408Q29| ; See also Ref. 11 

11. L.D. Paniak "Singular Matrix Models and the Kazakov-Migdal Model " M.Sc. 
Thesis, University of British Columbia, (in preparation). L.D. Paniak, N. Weiss, 
"Kazakov-Migdal Model with Logarithmic Potential and the Double Penner 
Matrix Model" ^lep-th 95010T7| 

12. S. Chauduri, H. Dykstra, J. Lykken, Mod. Phys. Lett A V6 1665 1991 



13 



